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1. Introduction 

In recent years, asymmetric simple exclusion processes (ASEP) have become a subject of 
increasing scientific interest because of their crucial role in the investigations of numerous 
dynamic phenomena in chemistry, physics and biology PHEj. They are important for 
understanding mechanisms of biopolymerization reptation polymer dynamics [3], 
diffusion through membrane channels [3], traffic problems jSj, and protein synthesis 
[31 El IE]- ASEP are non-equilibrium one-dimensional lattice models where particles 
interact only through hard-core exclusion and move preferentially in one direction. 

Although dynamic rules of asymmetric simple exclusion processes are very simple, 
they show a very rich dynamic phase behavior. One of the most striking features 
of ASEP with open boundaries is the occurrence of non-equilibrium phase transitions 
between stationary states that have no analogs in equilibrium systems El E] • The 
physics of these phase transitions can be explained by utilizing a phenomenological 
domain- wall theory According to this theory, the entrance, the bulk, and the 
exit of the system define their own stationary domains with specific uniform densities 
and currents. Domain walls exist in the border region between different domains, and 
the dynamics of these domain walls, which depends on the system parameters, in the 
stationary-state limit will determine the dominant phase. In the maximal-current phase, 
the steady-state density profile and current are due to the bulk dynamics, while the low- 
density (high-density) phase is enforced by the entrance (exit) rate. 

Recently, a new class of ASEP with particles occupying more than one lattice site 
has been investigated using various mean-field and continuum approaches |S] . These 
models provide a more realistic description of many biological processes such as RNA 
translation, vesicle locomotion along filaments, and proteins sliding along DNA. For 
example, in RNA translation a ribosome typically covers 10-12 codon sites but moves 
only one codon at a time [HI]- However, one more important dynamic feature, the non- 
uniformity of hopping rates, should also be considered for ASEP with extended objects 
in order to give a more realistic description of one- dimensional biological transport 
phenomena. This issue has been only briefly considered by using crude approximation 
methods and Monte Carlo simulations jSj. More thorough theoretical investigations 
of effects of inhomogeneity are needed. The main goal of this paper is to develop a 
theoretical and computational description of inhomogeneous ASEP with particles of 
arbitrary size. We consider here the totally asymmetric simple exclusion processes 
(TASEP), where particles can only move in one direction, although our theoretical 
arguments and our method can be extended to more general exclusion models. 

We present a theoretical study of the simplest inhomogeneous asymmetric exclusion 
model with particles of arbitrary size. In this model, the only nonuniform hopping rate 
is at the site which is in the middle of the lattice. We investigate this system by 
using a simple approximate theory, based on the domain-wall approach, and extensive 
computer simulations. The paper is organized as follows. In section 2 we outline our 
model and known results for homogeneous ASEP with extended objects, and we present 
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Figure 1. Model for TASEP with extended objects and a local inhomogeneity. System 
with size I = 3 particles is shown. Hopping rate in the bulk is 1. Hopping rate at special 
site k is q. Entrance and exit rates are a and f3, respectively. 

calculations from a simple approximate theory and domain- wall arguments. Monte 
Carlo simulations are discussed in section 3. We summarize and conclude in section 4. 

2. TASEP with extended objects 

2.1. Model 

We consider identical particles moving on a one- dimensional lattice consisting of iV sites, 
as shown in Fig. ^ Each particle may cover / > 1 sites. For convenience, we associate 
the position of each particle with the position of its left edge. In the bulk of the system, 
a particle located at site i will move to the next site i + 1 with the rate 1, given that 
site i + I is empty. Furthermore, when the particle is at the special site k, far away 
from the boundaries, the particle jumps forward to site k + 1 with the rate q (provided 
that site k + / is empty). Thus particles travel from the left to the right with uniform 
rates except at the special site k, which is the place of local inhomogeneity. Particles 
also can enter the system from the left with rate a if the first I sites are empty. When 
a particle reaches site N, it can exit with rate (3. We consider N to be very large, 
and in this thermodynamic limit the exact details of entrance and exit dynamics are 
not very important [7 . As in [S], we define the density of the particles in the system 
as the coverage density, namely, for M particles on the lattice the density is given by 
p = Ml/N . Since we discuss only the totally asymmetric simple exclusion model with 
extended objects, particles can move in only one direction, from left to right as shown 
in Fig. [T] 

2.2. Results for homogeneous TASEP with extended objects 

Analytical mean-field calculations and domain-wall arguments, supported by extensive 
Monte Carlo simulations, indicate that homogeneous TASEP with extended objects has 
the same three-phase diagram as homogeneous standard TASEP with I — 1 [7| [5] . 

For a > -jj^i and (3 > -/^j, the dynamics of the system is determined in the bulk, 
and we have a maximal- current phase with the stationary current and bulk density given 
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by 

i 1 ^ m 

J= (7TTTj^' p = 7TTi (1) 

When the particle entry is the rate-limiting step, which is realized for a < t=— j 
and a < j3, the system is in a low-density stationary phase with the following current 
and bulk densities: 

a(l — a) al 
l + (Z-l)a' P= l + {l-l)a ^' 

The conditions (3 < and /3 < a induce a high-density phase, which is governed 
by the exit dynamics. In this case the current and bulk densities are given by 

There are two types of phase transitions in this system. A first-order phase 
transition, with discontinuous change in density, is observed when a = (3. However, 
the transitions between low-density (high-density) and maximal-current phases are 
continuous. 



2.3. Approximate solutions of inhomogeneous TASEP with extended objects 

For simplicity, let us assume that the size of the system N is a large even number and 
that the special site with jumping rate q is positioned at k — N/2. The exact position 
of local inhomogeneity is not important as long as it is far away from both boundaries 

Hung. 

The special site k breaks the translational symmetry of original homogeneous 
TASEP. However, it also divides the system of size N into two homogeneous 
translationally invariant sublattices of size N/2. This observation allows us to consider 
our model with local inhomogeneity as two coupled homogeneous TASEP with extended 
objects. The sublattices are coupled by a condition that the stationary currents in the 
left and right subsystems should be the same. As a result, to calculate the properties 
of TASEP with extended objects and with local inhomogeneity, we can use the known 
results for homogeneous TASEP. This approach has been used successfully before in 
different inhomogeneous asymmetric exclusion models [TT| I12j. Note that the results of 
Ref. jTT] correspond to our special case 1 = 1, i.e., when each particle occupies only one 
site. 

The main assumption of our theoretical method is that TASEP with local 
inhomogeneity can be viewed as two independent homogeneous TASEP coupled only 
by the requirement to have the same steady-state currents. Then each of the 
sublattices may exist in one of three different stationary states, and there are nine 
possible phases in the overall system. However, since the particle current in the left 
subsystem should be equal to the current in the right subsystem, it is impossible to 
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have the maximal- current phase in one of the sublattices while the other sublattice 
is in the low-density or high-density phase. This observation eliminates four possible 
stationary phases and leaves only low-density/low-density (ld/ld), high-density/high- 
density (hd/hd), low-density/high-density (ld/hd), high-density/low-density (hd/ld) 
and maximal- current /maximal- current (mc/mc) stationary phases. In this notation, 
the first term corresponds to the state of the left sublattice, while the second term 
describes the right sublattice. 

Consider in a more detail the possibility of existence of an ld/hd phase in our 
system. This phase may exist only when a ^ (3 and 

a < — , < — . (4) 

Then, because the currents in the sublattices are the same, 

T _ q(i-q) (3(i-(3) (] 

l + (/-l)a l + (Z-l)/r {) 

which yields 

However, using the condition (0J) for a we may conclude that (3 > -yn~[-> which violates 
the condition for existence of the high-density phase on the right sublattice. Thus the 
ld/hd phase also cannot be realized in our system. 

The local inhomogeneity in our model is realized for different values of q. When 
q > 1 the fast jumping rate is introduced. In this case, the dynamics near the local 
inhomogeneity will not affect the overall dynamics in the system since crossing the 
special site k will not be a rate-limiting step. As a result, we have exactly the same 
phase diagram as for homogeneous TASEP with extended objects with ld/ld, hd/hd 
and mc/mc phases. The presence of local inhomogeneity will only modify the density 
profiles near the site k. 

The more interesting case is q < 1, when the dynamics near local inhomogeneity 
may determine the overall behavior of the system. Now consider in more detail the 
process of a particle crossing from the left sublattice to the right one. The particle first 
approaches the right sublattice when its left edge is on the site k — l + 1. Then it makes 
/ — 1 jumps with rate 1 and one jump with the rate q, if these motions are allowed, and 
the particle transfers completely into the right subsystem. Thus, the effective rate of 
moving from the left sublattice to the right sublattice is given by 

Qeff = l/ g +(Z-l) = l + (7) 
If the entrance to the left subsystem is the rate-limiting step, then the system will 
be in the ld/ld phase with the particle current and bulk density given in (J2J). This 
situation is realized for a less than some a*. Similarly, when the exit from the right 
subsystem determines the overall dynamics, the hd/hd phase exists with the particle 
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current and bulk density given in again for (3 < (3*. Parameters a* and (3* will be 
calculated explicitly below. 

When the transition from the left sublattice to the right sublattice becomes the rate- 
limiting step, the hd/ld phase is realized in our system. The parameters for existence 
of this phase can be found from the condition that the current through left sublattice 
is equal to the current through the right sublattice and equal to the current passing 
through the local inhomogeneity. Following the expressions (J2J) and ©, the stationary 
current and bulk density in the left sublattice are given by 

Peff(l-Peff) 



Jleft i + (3 eff (i-iy 

while in the right sublattice, 

a eff (l - a ef f) 



J, 



right 



l + a eff (l - 1)' 



Pleft 



Pright 



ft 



eff> 



la 



eff 



(8) 



(9) 



l + a eff (l- 1) 

where /3 e ff is the effective rate to exit from the left sublattice and a e /f is the effective 
rate to enter the right sublattice. The current passing through the local inhomogeneity 
can be written as 

Pleft (1 — Pright) 



Jo 



leff 



(10) 



I (1 — Pright + Pright /I) 

This expression can be understood as follows. The parameter q e ff, which is given 

explicitly in (j2J), is an effective rate of crossing the local inhomogeneity. The factor 

is the probability to find a particle at site k — I + 1, i.e., the particle is crossing from 



the left sublattice to the right sublattice. Finally, 



(1— pright) 



is the conditional 



(1 — Pright + Pright /I) 

probability for the corresponding sites in the right sublattice to be empty (as in jHj). 
The condition of stationary state implies that Ji e f t = J = J r i g ht, which yields 

_ Pleft ~ (1 - q e ff)[l - Pleft(l ~ 1)] 

Pright 



Qeff)[l~ Pleftil ~ 1)] 



Pright = Qeff Pleft- 

Solving these equations leads to expressions for the densities in the right and left 
sublattices: 



Pleft 



Pright 



(1 + Qeff)l 
2<7e//(J-l) 

Qeff Pie ft- 



1 



1 - 



\ '(1 + Vefff 



(12) 



The explicit formula for the current can be found from ()1()J) . This phase can only 
exist when a > /3 e ff and (3 > a e ff, and using (jHJ), © and (|T2*|) we obtain that hd/ld 
phase is stable for a > a* and (3 > (3*, where 



a 



f3* 



(1 + ge//)/ 

2? e //(i - 1) 



1 



V K^ + Qeff) 2 



(13) 



Our theoretical results can be easily checked in some special limiting cases. For 
q e ff = 1, and a > a*, (3 > (3*, we obtain that the bulk densities at the left and 
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right sublattices are the same and equal to pi e f t = p r i g ht = yft/(Vl + 1), i.e., the system 
is in the maximal- current phase, in agreement with known results [HE]. m the limit of 
q — > 0, the effective rate of crossing the local inhomogeneity is also very small, q e ff ~ Iq. 
In this limit the bulk densities reduce to pi e f t ~ 1 — q and p r i g ht ~ Z<Z, in agreement 
with intuitive expectations. These results can be understood as follows. Particles will 
move into the right sublattice with a rate of q, so that the spacing between the particles 
will be about l/q (distance traveled before the next particle comes into the sublattice), 
giving a coverage density of Iq. Now consider the left sublattice. It will be nearly full, so 
consider the backwards current of holes moving through the slow site. Holes will enter 
the left sublattice at a rate q, so the hole density will be q. Hence the coverage density 
in the left sublattice will be 1 — q. Another interesting limit is I 3> 1, which corresponds 
to the case of very large particles. In this limit we obtain J pa 1/7, as was found before 
for the homogeneous model [Zj. 

It is interesting to note that our theory predicts a "mixed" nature of phase 
transitions from hd/ld to ld/ld or hd/hd phases. At one sublattice the transition will be 
continuous, while at the other sublattice there is a jump in density, which corresponds 
to a first-order phase transition. 

3. Monte Carlo simulations and discussions 

Our theoretical approach gives correct results in limiting cases. However, to check 
the overall validity of our approximate theory we performed Monte Carlo computer 
simulations. 

Monte Carlo simulations were implemented similarly to those in jHj. Random 
sequential updating was used. Suppose that at some time the lattice contains M 
particles. Then in the next Monte Carlo step (MCS), M + 1 particles are chosen at 
random, in sequence, to attempt moves. They are selected from a pool containing the 
M particles on the lattice plus a free particle that can enter the system with probability 
a if the first I sites are empty. Particles on lattice site k advance with probability q 
if site k + I is empty. Particles on site N leave the system with probability (3. All 
other particles advance if they have room to move. Simulations were begun with an 
empty lattice and run until steady state was reached. For particles with I = 12, the size 
studied in greatest detail, the system size was N = 2500, and 3 x 10 5 MCS were allowed 
to reach steady state. For other particle sizes, the system was run for at least lOOiV 
MCS to reach steady state. After steady state was attained, systems were simulated for 
an additional 1.5 x 10 6 Monte Carlo steps, during which the current leaving the system 
was recorded constantly, and the particle density at each lattice site was sampled every 
100 MCS. Continuous time Monte Carlo [13] was used to reduce the computation time 
required. 

The phase diagram was determined in detail from simulations for I = 12, q = 0.2, 
and iV = 2500 (see Fig. |2J). From the theory, we expected the ld/ld phase to occur 
for a < a* = 0.1338, a < (3, the hd/hd phase to occur for (3 < (3* = 0.1338, (3 < a, 
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Figure 2. Phase diagram for I = 12, q = 0.2. Locations of phase transitions were 
determined via simulations to within 0.02 and were in agreement with predictions. 
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Figure 3. Sample density profiles from simulations for the ld/ld phase (top, a = 
0.05, j3 = 0.5), the hd/ld phase (middle, a = 0.5, j3 — 0.5), and the hd/hd phase 
(bottom, a = 0.5,(3 = 0.05). All cases have q = 0.2, I = 12, N = 2500. 



and the hd/ld phase to occur for a,/3 > 0.1338. Simulations were consistent with 
these expectations. Sample density profiles from each phase are shown in Fig. El The 
theoretical predictions for currents and bulk densities in each half of the system also 
matched the simulations to within 1% (see Fig 01) • 

The phase transitions were identified as follows. Near the transition between the 
hd/hd phase and the hd/ld phase, in the right half of the system, a e ff ~ (3. This 
situation is expected to produce shock waves in the right half of the system (see jHj). 
We detected these shock waves in simulations by observing the approximately linear 
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Figure 4. (a) Dependence of current J on a. (b) Dependence of bulk densities pi e ft 
(squares, solid curve) and p r i g ht (triangles, dashed curve) on a. Bulk densities were 
spatial and time averages over regions in which the density profile was approximately 
constant. pi e f t values are omitted near the phase transition because the density profile 
was approximately linear in the left half of the system. In all cases, symbols are 
simulation results, and curves are predicted values. Results are for j3 = 0.5, q = 0.2,1 = 
12, N = 2500. Error bars, determined from standard deviations for simulations run in 
triplicate, are smaller than size of symbols. 
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Figure 5. Sample density profile near the transition from the hd/hd to the hd/ld 
phase. Parameters are a = 0.2, q — 0.2,(3 — 0.13, Z = 12, N = 2500. The transition 
is expected to occur at (3 = 0.1338. Because of the proximity to the phase transition, 
the density profile is approximately linear in the right half of the system. 



density profile that they produce in the right half of the system. An example of a 
density profile for the hd/hd to hd/ld transition is shown in Fig. |SJ Similarly, the 
left half of the system exhibits an approximately linear density profile near the ld/ld 
to hd/ld transition, and both halves of the system have approximately linear density 
profiles near the ld/ld to hd/hd transition. Visual inspection of the density profiles 
allows the transition lines to be localized to within about 0.02 units. Other techniques 
would be required to find the transition lines more precisely. For example, we could 
determine the boundaries of the hd/ld phase by measuring the current as a function of 
a or (3 and determining at which a or (3 the current saturates. 

Simulations were conducted in each phase for other values of I (/ = 2, 5, 20) to 
confirm the theoretical predictions. System sizes N were always at least 200/, in keeping 
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with the assumption of large N. The theoretical predictions for the current and bulk 
densities were generally more accurate for larger particles. The greatest discrepancies 
between theory and simulations were seen for I = 2 in the hd/ld phase, but observed 
discrepancies were always less than 10%. These results show that our approximate 
theory becomes more exact for larger-size particles. This is because the particle-particle 
correlations decrease with increasing I. Consider a lattice consisting of iV = nl sites 
with M particles (each of size I). If we assume that particles are uniformly distributed, 
we can easily estimate the average number of empty sites x\ between two neighboring 
particles: x\ = nl ~^ Al = lx\. Thus the number of empty sites between the particles 
grows linearly with the particle size I. As result, particles correlate to lesser degree with 
each other, and our theoretical mean-field approach becomes more accurate. Similar 
conclusions have been reached in Ref. 0. 

The phase diagram and phase transitions in TASEP with extended objects and 
with a local inhomogeneity can also be easily explained with the help of the domain-wall 
theory j^j. As we discussed above, the introduction of local inhomogeneity divides the 
system into the two homogeneous sublattices, and the phase behavior in each sublattice 
is determined by the dynamics of domain walls. Thus there are always two domain 
walls present in the system. In the hd/ld phase at stationary state the domain walls 
are localized near the entrance and the exit of the system, while in ld/ld stationary 
phase the domain walls can be found diffusing near the local inhomogeneity and the 
exit. Similarly, in hd/hd stationary phase the domain walls are near the entrance and 
the local inhomogeneity. The phase transitions can be associated with the change of 
position of domain walls. For example, the transition from ld/ld phase to hd/ld phase 
(when q < 1) is the consequence of the motion of the domain wall from the position 
of local inhomogeneity to the entrance, while the second domain wall stays at the exit. 
This observation also explains the "mixed" nature of this phase transition. 

4. Summary and conclusions 

We investigated the effect of inhomogeneity in totally asymmetric simple exclusion 
processes with particles of arbitrary size. Specifically, the simple model with one 
inhomogeneous jumping rate in the middle of the system was considered. The model 
was solved using a simple approximate method. Our analytical approach was based 
on two assumptions. The first one is the fact that the local inhomogeneity divides the 
system into the two coupled homogeneous systems, for which very precise mean-field 
solutions already exist. In the second approximation we assumed that the dynamics of 
particles in each sublattice is independent, i.e., we neglected the correlations between 
the particles in the two coupled subsystems. Using these assumptions, we calculated 
explicitly the phase diagrams, bulk densities, and stationary currents. 

Our theoretical results indicate that for a fast jumping rate (q > 1), the phase 
diagrams, currents and bulk densities do not change as compared with homogeneous 
TASEP with extended objects. However, for the slow jumping rate, the situation is 
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different. As in the homogeneous TASEP, three stationary phases are found, although 
phase boundaries, stationary currents and bulk densities change significantly. Our 
theoretical predictions are well supported by extensive Monte Carlo simulations. Some 
small deviations between the analytical theory and computer simulations are due to the 
neglect of correlations near local inhomogeneity in our theoretical approach. However, 
the precision of our theoretical predictions increases with the size of the particles. The 
microscopic nature of phase diagram was analyzed using phenomenological domain- wall 
theory. 

Although we considered only totally asymmetric exclusion processes, our theoretical 
arguments can be easily extended to more general partially asymmetric exclusion 
processes where particles can hop in both directions. However, the important parameter 
of effective crossing rate from the left sublattice to the right should be calculated using 
more general expressions, g e // = V r > where r is mean first-passage time for the particle 
to move to the right sublattice. It can be calculated using standard procedures [TI] . 
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